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Summary 


Recent  activity  on  several  topics  in  propagation  through  atmospheric 
turbulence  is  described  in  this  report.  The  statistics  of  target  or  receiver 
irradiance,  conditioned  on  a prior  measurement  of  the  instantaneous  value, 
are  considered  under  conditions  of  moderately  strong  fluctuations,  and  it  is 
found  experimentally  that  such  scintillations  are  not  described  by  a bivariate 
log  normal  distribution.  Hence  as  turbulence  effects  become  stronger,  predic- 
tions of  the  state  of  the  turbulent  channel  at  a given  time  delay  following  a 
reference  measurement  will  be  subject  to  more  uncertainty.  However,  even 
though  the  exact  joint  distribution  is  not  analytically  known,  approximate 
predictions  based  on  the  bivariate  log  normal  are  quite  useful.  Such  predic- 
tions can  comprise  a basis  for  operating  a system  in  which  a low-power  laser 
beam  is  used  to  determine  a favorable  time  interval  for  firing  a high-power 
beam  in  turbulence. 

A.  study  of  computer  simulation  techniques  is  described  which  has  the 
gc  of  making  possible  realistic  models  of  propagation  paths  for  the  short- 
term, time-dependent  statistics  of  scintillation,  including  the  effects  of 
localized  or  intermittent  strength  of  turbulence.  It  is  concluded  that  the 
random  medium  must  be  more  completely  described  than  in  conventional  simula- 
tion approaches,  and  that  this  may  be  accomplished  through  an  empirical  trans- 
formation of  a gaussian  distribution  to  a distribution  which  realistically 
describes  the  highly  non-gaussian  behavior  of  refractive-index  fluctuations. 

An  expression  is  given  which  accomplishes  this  goal,  and  an  approach  for 
including  intermittency  is  described. 

The  important  future  topic  on  this  program  will  be  the  experimental  and 
analytical  determination  of  scintillations  from  laser-illuminated  diffuse 
and  general  targets.  A preliminary  analysis  of  scintillations  from  a diffuse 
target  is  given,  which  is  valid  in  the  near  field  of  the  target  illumination 
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I.  Introduction 


This  report  describes  recent  activity  on  several  topics  in  propagation 
through  atmospheric  turbulence.  In  Section  II  we  discuss  additional  aspects 
of  conditional  fading  statistics  to  supplement  the  treatment  in  an  earlier 
report.  In  Section  III  we  briefly  describe  work  on  the  useful  simulation 
of  short-term  scintillation  and  turbulence  intermittency  phenomena  on  a 
computer.  An  expanded  treatment  of  the  limiting  behavior  of  the  second 
moment  of  irradiance  is  given  in  Section  IV.  Finally,  in  Section  V we 
present  a preliminary  analysis  and  discussion  of  scintillations  from  diffuse 
and  general  sources,  which  will  comprise  the  primary  topic  for  future 
experimental  and  analytical  work  on  a follow-on  program. 

In  addition  to  the  work  reported  here,  we  have  continued  experiments  on 
target  irradiance  behavior  as  a function  of  wander-cancellation-tracking 
and  transmitter  and  turbulence  parameters.  In  particular,  this  work  has 
been  extended  to  10.6  microns.  Data  reduction  and  interpretation  is  being 
completed  and  a special  report  on  this  topic  will  be  issued.  We  are  also 
conducting  long-path,  strong-turbulence  scintillation  experiments  with  very 
small  receiver  apertures  and  wide  bandwidths;  this  work,  sponsored  by  the 
National  Science  Foundation,  builds  upon  earlier  efforts  on  the  present 
program  and  will  be  reported  later. 


II.  Further  Results  on  Conditional  Fading  Statistics 


In  the  preceding  report  on  this  program  , we  discussed  the  conditional 
statistics  of  scintillation  or  target  irradiance  fading  at  a time  t fol- 
lowing the  comparison  of  the  instantaneous  level  with  an  arbitrary  thresh- 
old. An  approach  was  given  which  is  readily  implemented  if  the  log  irra- 
diance fluctuations  have  a two-point  or  joint  normal  distribution,  and 
which  can  be  numerically  extended  to  a general  distribution. 

It  was  found  empirically  that  the  bivariate  normal  distribution  is 
accurate  for  weak  irradiance  fluctuations,  and  that  the  analytical  treat- 
ment of  the  conditional  statistics  is  successful  in  predicting  experimental 
results.  However,  it  was  mentioned  that  under  conditions  of  stronger 


1,  J,  R,  Kerr,  et  al,  "Propagation  of  Multiwavelenfcth  Laser  Radiation 
Through  Atmospheric  Turbulence",  RADC  Technical  Report,  April  1975. 
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fluctuations,  a significant  departure  from  bivariate  normality  occurs. 

The  purpose  of  this  section  is  to  elaborate  on  the  latter  point. 

We  confine  ourselves  to  strengths  of  fluctuations  below  t!»s  level  of 
"saturated  scintillations",2  which  will  be  treated  later.  It  will  in  fact 
be  shown  that  significant  departure  from  bivariate  normality  occurs  for 
moderate,  unsaturated  fluctuations. 

We  have  measured  the  conditional  distributions  of  the  log  irradiance 
fluctuations  from  a point  source  over  a range  of  turbulence  levels  (Table  1). 
For  each  experimental  run,  distributions  were  calculated  for  a number  of 
values  of  threshold  Uq)  and  time  interval  t,  using  400,000  data  points 
sampled  at  one  millisecond  intervals.  Autocorrelation  functions  were 
also  determined  as  in  Figure  6 of  Ref.  1. 


TABLE  1.  PARAMETERS  FOR  EXPERIMENTAL  DATA  ON  CONDITIONAL  FADING  DISTRIBUTIONS 


OVER  A 1 . 

6 km  PATH 

Run 

Mean  Wind  Speed 
(meters /sec) 

Log  Amplitude 

Variance  (0  2) 
x_z. 

Optical 

Wavelength  (p) 

1 

9 

0.0015 

10.6 

2 

3.5 

0.010 

10.6 

3 

9 

0.076 

0.488 

4 

1.5 

0.14 

10.6 

5 

1 

0.307 

10.6 

As  discussed  in  Ref.  1,  the  conditional  probability  distributions  for 

Run  ffl 

were  essentially  log  normal  for 

most  values  of 

threshold  and  time 

delay, 

and  the  observed  mean  and  variance  vs.  t for 

these  distributions 

agreed  with  the  analytical  predictions  based  on  this  assumption. 

Now  let  us  consider  a case  of  stronger  scintillation,  Run  #4.  Again 
the  conditional  distributions  appeared  qualitatively  (log)  normal,  as  sup 
ported  by  the  skewness  and  kurtosis  parameters  shown  in  Table  2. 


2*  1545’  Lawrence  and  J'  w-  Strohbehn,  Proc,  IEEE  58,  Oct.  1970,  1523- 
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fABLE  2.  * ^r™L®A^Li™XE%T)afis'lT0S  VARUSCE 

(-0,14). 


I 

— t (ms) 


-0.92 


0.075 


0.193 


-0.04 


0.013 


0.085 


0,055 


0.107 


0.048 


0.074 


0.048 


-2.6x10 


-6.6x10 


-3.6x10 


0.013 


However,  a comparison  of  observed  variances  and  those  T 

gausalan  assumption  (Table  3)  ahows  large  deviations;  1 , 

error  betweea  the  observed  and  predicted  normalized  variancea  la  132. 
error  betw  ,,  calculation  for  all  runs,  and 

in  Table  4 we  show  the  result  of  a similar  calcuiar 

ln  Tab  (o  2 > 0.1)  for  this  sampK 


1„  Table  4 we  show  the  resur.s  or  - ~ £or  this  sample 

it  is  seen  that  the  stronger- fluctuation  cases  ( x • ^ 

, significant  deviations  from  joint  normality.  It  s a so 
show  significant  ^ sensitive  to  departures 

predictions  such  as  those  in  Table  9 

fr°7 ;«:a"a„alysla  confirm,  the  non-ganaslan  behavior  of  the 

joint  attributions  In  the  stronger-fluctnatlon  runs.  This 
bat  in  the  weah-fluctuation  case  the  obaerved  variances  are  with  n the 
allowed  statistical  error  of  that  predicted  by  Joint  normality,  .bile 


- 1 , A.  >8,  .....  . i.„  Vv.... 


the*  strong-fluctuation  case,  the  deviations  are  outside  this  limit — 

especially  for  large  values  of  i . 

o 

A similar  analysis  of  saturated  log  irradiance  signals  will  be  carried 
out  in  the  future,  including  data  from  the  long-path  experiments  mentioned 
in  Section  I. 


TABLE  3.  COMPARISON  OF  PREDICTED  AND  OBSERVED  MEAN  OF  P(x(t  + r)/x(t)  « Jl 

a 

X 


— r(ms) 
X 

2 

4 

8 

16 

- .92 

Pred. 

- .856 

- .761 

- .544 

.215 

Obs. 

- .869 

- .77 

- .547 

.20 

i 

• 

o 

-p* 

Pred. 

.037 

.033 

.024 

.0094 

Obs. 

4x10”3 

.033 

.051 

.01 

1.07 

Pred. 

.995 

.885 

.632 

.250 

Obs. 

1.00 

.896 

.626 

.248 

2.G3 

Pred. 

1.89 

1.67 

1.19 

.475 

Obs. 

1.87 

1.54 

.962 

.314 

COMPARISON  OF 

PREDICTEL  AND  OBSERVED  VARIANCE 

OF  P(X(t  + r)/X(t) 

- i ) 

O 2 

0 

Jl 

X 

~ T (ms) 
X 

2 

4 

8 

16 

- .92 

Pred. 

.134 

.32 

.651 

.945 

Obs. 

.101 

.314 

.679 

.94 

i 

• 

o 

Pred. 

.034 

.32 

.651 

.945 

Obs. 

.088 

.255 

.615 

.96 

1.07 

Pred . 

.134 

.32 

.651 

.945 

Obs. 

.19 

.47 

.80  1 

.00 

2.03 

Pred. 

.134 

.32 

.651 

.945 

Obs. 

.268 

.603 

.95 

.99 

TABLE  4. 

RMS  ERROR  BETWEEN  OBSERVED  AND  PREDICTED  NORMALIZED  VARIANCES 

OF 

THE  CONDITIONAL 

LOG  IRRADIANCE  DISTRIBUTIONS 

Run 

a 2 
X 

, RMS  Error 

1 

0.0015 

0.04 

2 

0.010 

0.037 

3 

0.070 

0.027 

4 

0.14 

0.13 

5 

0.307 

0.11 
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III.  Computer  Simulation  for  Short-Term  Scintillation  Statistics  and 
Turbulence  Intermittency 

A.  Introduction 

The  goal  of  this  study  is  to  devise  computer  simulation  techniques 
which  will  make  it  possible  to  realistically  model  propagation  paths  for  the 
short-term,  time-dependent  statistics  of  scintillation,  including  the  effects 
of  localized  or  intermittent  strength  of  turbulence.  This  technique  comple- 
ments the  analytical  description  of  short-term  effects  given  in  preceding 
reports.1,3 

Obviously  the  most  ambitious  approach  is  to  model  the  complete 
refractive  index  field  over  the  path,  with  a differential  or  integral  formu- 
lation of  the  propagation  mechanism  which  yields  the  detailed  complex  field 
behavior  vs.  time  and  space.  Iteration  is  often  utilized  to  correctly 
include  the  multiple  scattering  regime.  There  are  two  primary  difficulties 
with  this  approach: 

1.,  Computation  requirements  are  prohibitive. 

2.  The  refractive  field,  having  nongaussian  statistics, 

is  not  readily  modelled  in  detail. 

Let  us  consider  concrete  examples.  Over  a limited  pathlength, 

perturbation  analysis  yields  the  following  form  for  the  complex  phase  with 

4 

an  elementary  source: 


4>(x,y,L)  = 


U, (dK  ,dK  ,L)= 
<J>  x y 


g ;x-x\y-y'  .L-z'^Cx'  ,y'  ,z') 


dz'  G(K.L-z')  U (dK  ,dK  ,z') 

'v*  e x y 


(1) 


where  e is  the 
1 


(stochastic) 


dielectric  perturbation,  and  the  second 


3.  J.  R.  Kerr,  et  al,  "Propagation  of  Multiwavelength  Laser  Radiation 
Through  Atmospheric  Turbulence",  RADC  Technical  Report,  November  1974. 

4.  V.I.  Tatarskii,  The  Effects  of  the  Turbulent  Atmosphere  on  Wave  Propaga- 
tion,  available  from  National  Technical  Information  Service  (#TT-68-50464) 
1971. 
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expression  is  simply  the  Fourier-Stieljes  form  of  the  first.  In  principle, 
could  be  modelled  in  detail  in  three  dimensions  and  with  time  evolution; 
the  statistics  of  * would  then  b . determined  from  a number  of  such  realiza- 
tions; „ i.e.,  a limited  ensemble.  Analytically,  the  solution  is  usually  con- 
fined to  a statistical  quantity  such  as  an  autocorrelation  function,  which 
involves  partial  Fourier  transforms  and  Kolmogorov  spectra,  so  that  the 
detailed  knowledge  of  the  dielectric  or  optical  field  for  a particular  case 
is  not  involved. 

B.  Use  of  Short-Path  Phase  Function 

To  model  the  system  represented  by  Eq.  (1),  we  must  generate 

realizations  of  the  three-dimensional  dielectric  field  e. (x' ,y' ,z') . How- 

3 1 

ever,  this  random  variable  is  markedly  nongaussian,  and  therefore  unsuitable 
for  conventional  simulation  techniques.  The  approach  generally  used"*’*’  is  to 
derive  a much  simpler  expression  like  Eq.  (1),  for  the  real  phase  only,  as 
follows.  Let  U be  the  complex  amplitude,  and  write 

U(x,y,z)  - er(x’y‘2)w(x,y,z)  (2) 

where  T is  the  real  phase  and  is  assumed  to  have  the  short-path  solution 
z o + Az 

r - I dz  e1(x,y,z)  • (3) 

z 

o 

When  U from  Eq.  (2)  is  put  back  into  the  wave  equation,  it  is  found  that 
w satisfies 

I” i2k  — + e~F  72  eF  1 w = 0 . (4) 

|_  Bz  xy  J 

It  is  then  assumed  that  a sufficient  number  of  independent  regions  are 
contained  in  the  integral  (3)  to  render  T a gaussian  random  variable,  and 

5.  J.  Herrman  and  J.  C.  Bradley,  Numerical  Calculation  of  Light  Propagation, 
Lincoln  Laboratory  Rept.  LTP-10,  July  1971. 

6.  W.  P.  Brown,  Jr.,  High  Energy  Laser  Propagation,  Midterm  Technical  Report, 
Contract  N00014-73-C-0460,  Hughes  Research  Laboratories,  September  1973. 
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the  modelling  is  limited  to  r(x,y)  realizations  with  independent  (uncor- 
related) iterations  for  each  pathlength  increment  (Az) . Although  Eq.  (3) 
is  valid  only  over  very  limited  path’-ngths,  the  "return  to  the  differential 
equation"  (4)  at  every  Az  iteration  l*ads  to  results  which  are  valid  for 
arbitrary  pathlengths  nd  which  therefore  include  the  saturation  of  scintil- 
2 

lations. 

Finally,  the  generation  of  particular  r<x,y,zo>  realizations  with  the 
required  transverse  correlation  function  is  accomplished  through  the  follow- 
ing principle:  a multidimensional  gaussian  random  variable  f(jr)  may  be 

generated  by  convolving  a gaussian  white  noise  function  h(^)  with  a ieter- 
ministic  function  e(^)  having  the  desired  autocorrelation  for  f : 

£<*>  - «<JE>  * <5> 

In  practice,  one  usually  employs  the  spectral  form  of  this  expression: 


/E(jjC)  H(JC)  e1^  £ 


(6) 


where  E(£)  is  the  spatial  power  spectrum  desired  for  f(^).  In  practice, 
the  problem  is  discretized,  so  that  H represents^  grid  of  gaussian  random 
numbers;  E(£)  is  the  familiar  Kolmogorov  spectrum  and  Hx^)  has  the 
usual  transverse  phase  structure  function  for  propagation  through  turbulence. 

Finally,  the  above  approach  may  be  extended  to  time  dependent  behavior 
by  invoking  the  frozen-in  turbulence  or  Taylor  hypothesis  ’ and  having  the 
randomly-generated  grid  move  with  the  postulated  wind  vector. 


C.  Possible  Extensions  of  the  Above  Technique 

The  detailed  nature  of  the  ^ field  is  dropped  out  of  the  above 
technique  by  integrating  over  short,  discrete  paths  and  assuming  that  the 
result  is  gaussian  and  uncorrelated  in  adlacent  z-segments,  even  though  the 
integrand  is  nongaussian.  For  Az  much  longer  than  the  outer  scale  of 
turbulence2’4  but  short  enough  for  Eq.  (2)  to  be  valid,  this  should  certainly 
give  accurate  statistical  results  over  a reasonable  ensemble  of  U realizations 


In  our  study  of  short-term  statistics,  including  turbulence 
intermittency,  we  have  been  interested  in  the  detailed  statistics  of  the 
turbulence  field  and  in  the  detailed  time  behavior  of  the  scintilla i ig 
irradiance,  including  joint  distributions  and  conditional  statistics  as  in 
Section  II.  Such  details  are  obviously  lost  by  integrating  over  pathlength 

segments  before  generating  a realization. 

In  order  to  extend  the  above  or  any  other  coding  technique  to 
include  the  detailed  e ^ statistics,  we  suggest  the  use  of  a versatile 
"Johnson"  trensformation7  which  we  have  found  can  adequately  represent  the 
nongaussian  variable  in  terms  of  a gaussian  variable  r^: 


ni 

ei  * ci sinh  r2 


(7) 


where  and  C are  determined  by  the  first  two  momemts  of  This 

technique  simply  comprises  an  empirical  fit  of  the  probability  distribution 
of  a nongaussian  process  to  the  gaussian  distribution,  obviously  with  the 
loss  of  independent  information  on  higher  moments.  The  autocorrelation 
or  spectrum  (E)  of  must  be  calculated  from  the  known  (isotropic)  form 

for  e^. 

We  thus  write  the  generating  expression  for  T,  from  Eqs.  (3), 
(5),  and  (7),  as 


T(x 


ZQ  + Az  i “ 

,y)  = ^ | dz  C±  sinh  f j jj  dx'dy'dz'  h(x-x' 


-x' ,y-y' , z-z ' ) 


e(x' ,y ' ,z') 


(8) 


The  segments  Az  should  now  be  much  less  than  the  outer  scale,  in  order 
to  preserve  maximum  detail.  This  is  also  true  of  the  transverse  discretizing 
or  grid  scale.  The  transformed  refractive  field  may  again  be  translated 


7.  G.  F.  Hahn  and  S.  F.  Shapiro,  Statistical  Models  in  Engineering,  Wilev, 
New  York,  ]967. 


with  the  wind,  in  order  to  determine  time  behavior. 

In  principle,  this  technique  may  permit  the  determination  of  the 
detailed  statistics  of  scintillation,  including  two-point  or  conditional 
behavior,  higher  moments  (degree  of  fit  to  a log  normal),  autocorrelation 
functions,  integral  scales,  and  averaging  time  effects  on  spread  in  measured 
quantities  such  as  log  amplitude  variance.  However,  it  is  clear  that  the 
required  computer  time  may  be  prohibitive,  This  is  currently  under  inves- 
tigation for  a two-dimensional  representation  (x,z), 

A second  extension  of  the  technique  described  in  Sec,  III-B  is 
to  postulate  large-scale  random  variations  of  turbulence  strength  (envelope 
of  the  fluctuations),  representing  intermittency  of  turbulence.  The 

nature  of  such  variations  has  been  discussed  in  Ref,  3,  The  Johnson  trans- 
formation and  attendant  detail  is  not  required  t<  r the  basic  results  of 
interest,  and  in  fact,  as  discussed  in  Ref,  3,  the  intermittency  is  only 
manifested  in  higher  moments  of  which  are  lost  in  the  empirical  trans- 
formation process, 

IV.  Asymptotic  Expansion  of  the  Second  Moment  of  Irradiance 

Using  physical  reasoning,  we  have  argued  in  recent  reports  that  the 

normalized  variance  of  irradiance  (c^2)  should  approach  unity  at  large  values 

of  transmitter  diameter  (D)  over  coherence  length  (p  ) , where  p is  deter- 

o o 

mined  over  the  reciprocal  (target-to-transmitter)  path  through  turbulence. 

This  was  shown  directly  from  the  analytic  expression  for  the  fourth  moment 
of  amplitude  in  Ref.  3. 

Efforts  have  continued  to  numerically  or  analytically  derive  further 
details  of  the  curve  of  cr^.2  vs.  (D/p^),  in  order  to  complement  the  phenomeno- 
logical treatment  which  successfully  predicts  the  essential  features.®  The 

9 

same  problem  has  been  treated  by  Banakh,  et  al,  using  certain  approximations 
and  a Monte  Carlo  technique,  with  results  which  are  highly  useful  but  which 
have  certain  limitations  as  discussed  in  Ref.  8. 

The  purpose  of  this  section  is  to  extend  the  asymptotic  analysis  of 
Ref.  3,  using  an  inverse  expansion  of  (D/pq).  We  begin  with  the  integral 

8.  J.  R.  Kerr,  et  al,  "Propagation  of  Multiwavelength  Laser  Radiation  Through 
Atmospheric  Turbulence",  RADC  Technical  Report,  May  1974. 

9.  V.  A.  Banakh,  et  al,  J.  Opt.  Soc.  Am.  64^  516-518,  April  1974. 


11- 


expression  for  the  fourth  amplitude  moment  (second  moment  of  irradiance) : 


2 k 

< — 

(2ttL)/ 


llll  dhd^: 


3d£^  exp 


2 , 2 . 2 2 
P1  + p2  + p3  + p4 

D2/2 


• exp 


/ n \ 5/3  T 5/3  5/3 

(p^)  ^73  j^l"^  + l£3_£J 


5/3  5/3  5/3 

+ + |£2-p3I  - i£2-£*l  ■ |£ 


r£2 1 Ji 


(9) 


We  let  x^  ■ £i/d  and  ® “ (D^Po)  and  rewrite  (9) 


as 


< j2>  " ( ml)  HU  d*id*2d*3d*4  «P  {-  *<x*  + *2  + x3  + } 


[ ^ 


5/3 


5/3 


exp  <j-  B fag-22 1 ♦ li3“*4 1 + I*!-: 


5/3 


*4 1 


5/3  5/3 

+ 1 2L2  — — 3 1 - I *2-*J 


5/3] 

- lij-Sjl  } 


(10) 


For  large  B the  integrand  is  significantly  different  from  zero  only 
near  the  points  x^  “ 2L2  and  = Consequently,  one  should  be  able  to 

obtain  approximations  to  the  integral  by  considering  contributions  from  near 
these  two  points.  We  note  further  that  the  integrand  is  symmetric  with 
respect  to  interchange  of  Xj  and  x^  which  means  that  we  need  only  to  calcu- 
late  the  contributions  to  the  integral  from  x-^  “ x.2  and  multiply  this  result 
by  two. 

The  following  changes  of  variable  simplify  the  calculation.  Let 
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so  that  x^-x^  = r+s+t;  x^x^  " s+t  * and  x^-x.3  = These  substitutions 

transform  the  integral  to 


(JJ| 


dx2dt^dsdr  exp 


2 2 
8x2  + 8x2t  cos0t  - 4t 


- 2r^  - 2s^  - 4x„r  cos0  + 4x0s  cos0  - 4st  cos(0  -0  )i 
2 r 2 s 8 t J 

• exp  | - 6rS/3  -es5/3  - 6t5/3  - 6|r+4+t|5/3  + 6|e«|5/3 

+ 6 J s+t | 5/3 1 (12) 

where  x2  has  been  taken  to  be  the  polar  axis  and  0r,  0g,  and  0fc  are  the 

angles  between  x.2  and  _r,£,_t  respectively. 

We  note  that  for  8>>1  the  major  contribution  to  the  integral  comes  from 

near  the  points  r * 0,  s = 0.  This  corresponds  to  the  points  x^  * x.2» 

Xg  ■ *The  factor  exp  | 8t"*^  - 8|r+s+t|^^  + 3|r+t|^^  + 8|s+t|^^  | is 

bounded  and  does  not  extend  the  effective  range  of  the  r and  s integra- 

-3/5 

tions.  This  effective  range  is  0(8  ).  On  the  other  hand  the  t inte- 

gration is  not  so  restricted  since  the  above  factor  approaches  unity  as 

t -►  “ for  finite  r and  s.  The  range  of  the  t integration  is  0(1)  due  to 
-4t2 

the  factor  e 

We  will  now  exploit  the  narrow  range  of  the  rand  s integrations  by 

_6r5/3  j -gs5/3 

expanding  the  integrand, apart  from  the  factors  re  and  se  , in  a 

Taylor  series  about  r = s = 0,  Since 
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J drr 


5r  V — > 6 


This  will  result  in  an  expansion  in  inverse  powers  of  6.  We  thus  write 


< I2> 


(!!l)  jjd*2d^  6XP  {"  8X2  + ®X2t  CO80t  ' 4t2}  Q(V-^ 


where 


Q(x2»t.)  = drd  ij  exp(-0r5^3~Bs5^3)  exp  j -2r2-2s2-4x2r  cost>r 

+ 4x„s  cos0  - 4st  eos(0  -0.)>  exp  (-gf(r,s,t))  (1 

2 s s t ) 


f (r,s,t)=  t5/3  + i iHs+t  | 5/3  - |r+tj5/3  - |s+jt  | 5/3 


Now  expanding  f(r,s,t)  in  a Taylor's  series  about  r - s = 0,  we  have 


f (r,s,t)=  f (0,0,t)  + r 


, 2 2. 

8f  8f  . r_  jT£ 

8r  - 8 8s  . 2!  a 2 _ 

1 r=s=0  r=s=0  8r  1 r=s=0 


2 2 2 
2rs  8Zf  s 8 f 

2!  8r8s  n s!  . 2|  . 

r*s=0  8s  ' r l>  0 


We  then  note  the  following  equalities: 
f (0,s,t)=  f(r,0,t)  = 0 


;'***  V*->  • -'V  '%£  * f 


0 


3f 

9r 


3f  I 


r=s=0 


9s 


r=s=0 


92f 


9r 


92f 


r=s=0  9s  1 r=s=0 


2| 


92f 


9r9s 


r=s=0 


J t ^3  (cos(0r-0s)  - J cos (0^-0 cos (0^-0^)  ^ 


93f 


9r' 


93f 


r=s=0 


9s' 


r=s=G 


93f 


9r29s 


r=s*0 


~ t ^3  cos(0  -0  )cos(0  -0  ) - t cos (0  -0  ) 
9 rt  r s 9 st 


35  -4/3 


27 


cos  (0  -0  )cos(0  -0  ) 
r t st 


and 


93f 


9r92s 


r=s=0 


t ^3  cos(0g-0t)cos(0r-0g)  - |-  t"^3  cos(0r~0t) 


35  -4/3 


' 27 


cos  (0  -0  )cos(0  -0  ) 
sc  r t 


(17) 


Therefore 


exp(-0f (r ,s,t))  =1-4  8rst  ”^3  (cos(0  -0  ) - 4 cos(0  -0  ) cos(0  -0 

J v rs3  rt  st/ 


— 0r2st  ^3  cos(0  -0  ) cos(0  -0  ) - 4 cos (0  -0  ) 

L 9 r t r s 9 s t 


3! 


+ 44  cos2(0  -0  ) cos(0  -0  )1 
27  r t s t J 

3 2 -4/3  T 10  <; 

-■jjgrs  t [— g-  cos(0g-0t)  cos(0r-0s)  - ~ cos(0r~0t) 


35 


27  cos  (0g-0t)  cos(0r-0t)J  - 
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+ 1 25  „2_2_2^-2/3 


cos 


2!  9 e2r2®2t  2 3 [c°s(er-eg) 

<vv]2 

+ il(f ) B4rVt"4/3  [coS(8r-es) 

<VV]4 


- cos (6r~ e t ) 


jcos(er-et) 


cos 


+ i!(|)e6r6s6t"2  [cos(er-es) 


j cos(6r-0t) 


cos 


<es“et)]< 


(18) 


Although  (18)  does  not  'include  all  teruit  up  to  O(r^s^),  those  neglected 

either  contribute  zero  under  the  d0^  and  d&g  integrations  (see  below)  or 

result  in  terms  of  lower  order  in  8.  The  last  term  listed  is  proportional 
-2 

to  t and  hence  the  t integration  of  this  term  diverges.  This  is  because 

our  expansion  emphasizes  the  region  near  t = 0.  The  contribution  of  this 

—18 /5 

term  turns  out  to  be  0(8  ) so  that  there  is  ,.o  point  in  including  higher 

order'  terms. 

As  mentioned  above,  the  d0^  amd  d0g  integrations  drive  many  terms  in  the 

expansion  to  zero,  namely  all  those  which  are  periodic  with  period  2tt. 

-14/5 

Keeping  only  those  terms  which  are  non-zero  and  of  0(8  ) and  greater: 


Q(x2»I> 


2tt  2tt  « « 

/ d6rJ  d0„J  drrJ 


dss  exp(-8r 


-Br5'3-*,5'3 


• { 1 + 7 (f ) B2rVt'2/3  [cos(er-es)  - i cos(8t-0c) 

• cos(0g-©t)j2 
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+ \\  (f)  ' sVA-4'3  [coS<er-e3)  - i cos(er-8,) 


cos(8s-8t)]4  } 


+ 0(8-18/5) 


The  factor  2 appears  since  we  are  considering  only  contributions  to  < I > 

from  near  x.^  = — 2 ant^  not  ^rom  = 25^  (—  £ + Jl  = ^ • 

Performing  the  angle  integrations  gives 


J»00  m <30 

drr  I dss  exp(-9r5^3-Bs5^3)  | 4ir2  + 

o *0 


2 + 13/ 
9 


1 /5\2  „2  2 2 -2/3  121  2 1 / 5 \ 4 .4  4 4 -4/3  \ 

2 (3)  6 r S 1 + 144  * ‘ 4!  (3)  B T S t j 


+ 0(B-18/5) 


3,2.  (3)  2 r2(i)6-12/5  + l^  r^)  6-14/5t-2/3 

+ i|i_8  ,2  (5}2  r2(^}  6-16/5t-A/3  + 0(6-18/5) 


Since  to  this  order  Q(x2,t)  is  independent  of  x_2  and  d6t>  we  may  make  the 


substitution  v = 2x2~t  so  that 


2 2 2 

v = 4x2  - 4x2tcos9t  + t 

and  dx.2  = ^ dv  = -^  vdvd9v 


resulting  in 


<I2  > 


(8l)  \ | d0t  { dttJ  dev  J dv  v exp(-2v2-2t2) 


J J dt  t exp(-2t2)  | 8tt2  (|)  r2(|)  g' 
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00 
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T 
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wt 

( 

3 
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00 

r 

/W)  X 
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< 

| dt 

3 

13ir2 
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} + 0(B-1S/5 

(^)4{ 

2 V| 
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r2(f)  3-/3  r( 


V\f)  8-14/5  + 


-1/3 


12  \ -14/5  , 121  2 _4  Ul\  ( 5 


1728  8 
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b-16/5}+  0((f18/5) 


The  mean  irradiance  Is 


I > 


’ (Sl)  j j d£idfi.; 


exp 


2 2 
P1  + p2 


2 

D /2 


exp 


<£) 


5/3 


0J|  NJ 


’ n5/3 


5/3 


(sl)  J[%d^2  <“P(-2Xj-2Xj-e|x1-X2 


1 5/3 


exp(-4x2-4x2r  cos6-2r2)exp(-Br5^3)  (23) 


As  above,  we  let  v •»  2x„  + r_,  giving 


< I > 


d_r  exp(-r2-(3r3^3) 


2 2 C 

i I 

\ 2ttL  / 4 J 


(Hl)  T { dr  r exp(-Br5/3)  j 1-r2  + ~ - . . . j 


(§e)  r f r(t)  e'6/5  - ^12/5> 


(24) 


The  second  central  moment  of  irradiance  normalized  to  the  mean  irra— 
d lance  is  then 


r2-^f  C rd)  (f) 


• r'2(l)  B'2/5  + lfk  2-4/3  r(j)  (})*  r2(f)  r-2(f) 


• B-4/5  + 0(e-6/5) 


1 + 3.132B_2/5  + 9.414g_4/5  + 0(B_6/5) 


(25) 
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We  may  thus  write 


(1)  al 

= 1 + 3.1323”2'5 

(26a) 

(2)  al 

- 1 + 3.1323_2/5  + 9.4143_4/5 

(26b) 

2 

as  the  one-  and  two-term  asymptotic  approximation  to  . The  one-term 
approximation  agrees  exactly  with  the  result  obtained  by  a different  method 
by  Gochelashvily.^ 

(1)  0^2  an<j  (2) 0 * are  tabulated  for  a range  of  values  of  P in 

Table  5,  along  with  D/po  and  Dg  used  by  the  Russian  investigators.  Curves 
are  shown  vs.  /D^  in  Figure  1,  along  with  the  result  from  Ref.  S.  We  note 
that  our  result,  which  should  be  accurate  for  large  3 (Dg) » does  not  abree 
well  with  the  Monte-Carlo  (M-C)  integration  f /t  /Dg  > 20.  We  surmise  there- 
fore that  our  results  represent  a useful  modit.  .ation  at  large  Dg. 

The  M-C  calculation  should  be  more  accurate  at  small  Dg  or  3,  and  it 
agrees  with  the  one-term  expansion  for  8 < ^ < 20.  However,  the  two- 

term  expansion  appears  to  fit  the  data  somewhat  better.  ^ 

Two  final  points  can  be  made.  First,  the  structure  function  (^r  ) is 

valid  only  for  path  lengths  (L)  such  that  t'c<<L«Li»  where  is  the  propa- 
gation distance  at  which  the  near  field  decays  to  e of  its  value  at  the 
aperture  and  h±  is  the  propagation  distance  at  which  the  coherence  length 
js  of  the  order  of  the  inner  scale. ^ For  example,  if  'v  pQ  'v  1 mm  and 
D v 10  cm,  we  have  D/pq  = 100,  3 = 2150,  and  ^ = 66.  fhe  above  analysis 
is  then  valid  for  3«2150,  while  for  3»2150  the  structure  function  is  pro- 
portional to  r2  and  the  resulting  new  integral  can  be  evaluated  by  the  same 
method  used  above. 


10.  Prokhorov,  A.M.,  et  al,  Proc.  IEEE  63,  790-811,  May  1975. 

11.  Lutomirski,  R.F.  and  Yura,  H.T. , J.  Opt.  Soc.  Am.  61,  482-487,  April 
1971. 
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TABLE  5. 

TABULATION 

OF  NORMALIZED 
TERM 

VARIANCE  OF  IRRADIANCE  FOR  ONE- 
EXPANSIONS 

AND  TWO' 

D/oo 

8 

f2'5 

(Doj2 

(2)0l 

4 

3.48 

8 

0.435 

2.363 

4.147 

6 

5.66 

18 

0.3147 

1.986 

2.918 

8 

8.00 

32 

0.250 

1.783 

2.3/1 

10 

10.45 

50 

0.2091 

1.655 

2.067 

12 

13.01 

72 

0.1807 

1.566 

1.874 

14 

15.66 

98 

0.1598 

1.500 

1.741 

16 

18.38 

128 

0.1436 

1.450 

1.644 

18 

21.17 

162 

0.1307 

1.409 

1.570 

20 

24.02 

200 

0.1201 

1.376 

1.512 

22 

26.93 

242 

0.1113 

1.349 

1.465 

24 

29.90 

288 

0.1038 

1.325 

1.427 

26 

32.91 

338 

0.09737 

1.305 

1 394 

28 

35.97 

392 

0.09177 

1.287 

1.367 

30 

39.08 

450 

0.08684 

1.272 

1.343 

32 

42.22 

512 

0.08247 

1.258 

1.322 

34 

45.41 

578 

0.07856 

1.246 

1.304 

36 

48.63 

648 

0.07505 

1.235 

1.288 

38 

51.89 

722 

0.07188 

1.225 

1.274 

40 

55.19 

800 

0.06899 

1.216 

1.261 

780 

1949.4 

304200 

0.00641 

1.020 

1.020 

Note:  fi='(D/po)  ' - y Dg 

9 2 S/1 

D =1.095C  \ L(D)J/J 

S T1 
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The  second  point  is  that  the  reason  our  asymptotic  series^breaks^down 
after  the  second  term  is  that  expanding  quantities  like  |r+t|  as  t times 

a series  in  powers  of  r/t  is  valid  only  if  r«t.  The  effective  range  of 

the  r integration  is  0(g"3/5)  which  means  that  if  t»B  , the  expansions 

are  valid.  Presumably  we  could  avoid  the  divergence  caused  by  the  factors 
of  - as  t -*•  0 ty  cutting  off  the  integration  at  e.g.  t=  3 • This  intro- 

duces an  error  of  0(ff6/5),  so  that  terms  up  to  ff4/3  should  be  acceptable. 

The  t = 0 region  corresponds  to  the  region  where  x£  * x^  i.e.,  where  the 
two  maxima  of  the  integrand  coincide.  It  is  not  possible  to  handle  the 
integration  accurately  in  this  case,  and  as  6 gets  smaller  this  error 
becomes  more  significant. 

V.  Scintillations  from  Diffuse  and  General  Sources— A Preliminary  Analysis 
and  Discussion 

The  major  new  area  of  investigation  on  this  program  is  that  of  turbu- 
lence effects  on  radiation  from  diffuse  and  partially  coherent  sources,  with 
particular  attention  to  laser-illuminated  targets  and  the  performance  of 
coherent  optical  adaptive  (COAT)  systems.  This  represents  a natural  exten- 
sion of  the  recent  work  done  on  focused  wander-tracking  transmitters.  The 
experimental  effort  will  involve  the  determination  of  scintillation  statistics 
from  diffuse  targets  (speckles)  and  glints,  including  infrared  wavelengths  of 

primary  interest. 

The  purpose  of  this  section  is  to  present  a preliminary  analysis  for 
the  case  of  a diffuse  target,  coherently  illuminated,  with  the  point  of 
observation  in  the  near  field  of  the  target  illumination.  This  treatment 
includes  the  effects  of  turbulence  between  the  illuminator  and  target  and 
on  the  return  path  to  the  observer.  Generalizations  of  the  illuminator  and 
target  characteristics  and  geometry  are  underway  and  will  be  reported  in  the 

future. 

A.  Mean  Irradiance  at  Receiver 

The  source,  target,  and  receiver  configuration  is  shown  in  Figure  2. 
The  present  analysis  is  confined  to  the  case  of  a TEMQ{)  laser  illuminator. 

The  source  and  target  are  assumed  to  be  much  smaller  than  the  path  length 
(L),  and  the  distance  between  the  receiver  and  source  is  greater  than  the 
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Figure  2.  Target-Illuminator  Configuration 
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source  size  and  much  smaller  than  the  path  length.  These  geometric  condi- 
tions confine  the  problem  to  small  angles  and  ensures  that  the  outgoing  and 
returning  radiation  experience  independent  turbulence  regions;  the  latter 
limitation  is  thought  to  be  inessential  owing  to  the  diffuse  target  charac- 
teristics and  will  be  relaxed  in  subsequent  work. 

We  write  the  source  amplitude  distribution  as 


U (r)“  U exp 
o — o 


ikr2  \ 
2F  ) 


(27) 


where  a and  F are  the  characteristics  beam  radius  and  focal  length, 
o 

respectively.  The  field  at  the  target  is  written  from  the  extended  Huygens- 

, , J , 12,13 

Fresnel  principle  as 


U(p) 


, ikL 
ke 

2iriL 


f U (r)  exp  I"  lk  + M£.,r)  dr 

Jo-  2L  1 J 


(28) 


where  4^  describes  the  effects  of  the  inhomogeneous  medium  on  the  propaga- 
tion of  a spherical  wave.  Combining  (27)  and  (28),  we  have 


U(£) 
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kelk  [L  + 2L  ]«o 

fexp  T-  ^ + — 

(l  ,i),2 

2*iL 

J p 1 2« 2 2L 

V T ) 
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+ 4^  (£,r)  dr 


ik 

L 


£’£ 


(29) 


In  particular,  this  applies  to  the  special  cases  of  a focused  (L=F)  colli- 
mated (F  -*■  “)  beam  respectivel. 

The  field  at  the  receiver  is  written  by  reapplying  the  Huygens-Fresnel 
principle  to  the  field  at  the  target: 


12.  Lutomirski,  R.F.  and  Yura,  II. T.,  Applied  Optics  10,  1652-1658,  July  1971. 

13.  Yura,  H. , Applied  Optics  JJ,  1399-1406,  June  1972, 
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ke 


U(o> 


ikf"L  + £2.1 

L 2L J 

2rriL 


J U'(p)  exp  (p2-2£*£)  + i|*2(£,£)J  d£ 


(30) 


where  U'(p_)  is  the  field  solution  after  the  reflection  from  the  target, 
and  \f>  represents  the  turbulence  effect  from  the  target  to  the  receiver. 

The  mean  intensity  at  the  receiver  is  then 


2 

<I(£)>  = <|u(£)|2>  = (2^)  Jj  dp^  <u'(p1)u'*(p2)> 

* exy  [f  ((p?  ~ p22)  - 2£*<pr£2))] 

< exp  r>)J2(£»£1)  + ^2*(£*£2)J  > (31) 


Through  the  assumption  of  a diffuse  target,  the  reflected  beam  suffers 
a random  phase  delay  from  point-to-point  over  the  target,  and  we  write 


<U,(Ll)U'*(p2)>  = <I(PX)>  <5 (P^-P,,) 


(32) 


Using  this  in  Eq.  (31),  the  mean  intensity  becomes 


d(£)>  - (2^;)  J d£x  <|u(p1)|2>  <exp[iP2(£,£1)  + '^2*(£,£1)]: 


(33) 


where  the  mean  exponential  term  is  unity  from  considerations  of  energy  con- 
13 

servation. 


-26- 


The  resultant  mean  Intensity  at  the  receiver  is  then  simply 


(34) 


To  complete  the  solution,  we  use  Eq,  (29)  with  Eq.  (34).  We  note 

2 

that  the  structure  function  gives  us 


<exp  [^(p.rj)  + 'l'1*(P,r2)J>  = e 


<t) 


5/3 


(35) 


For  the  focused  beam,  we  then  have  (r  = |r.  -r  |): 


:|U(p)|' 


(m;)2  °o2  JJ  dhdI-2 


exp 


2.2 

rl+r2 

2a2 


11c 
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(— 1~— 2} 


(36) 


Carrying  out  the  integration  indicated  in  (34),  involving  the  Fourier-Bessel 
integral,  we  have  finally 


<!(£)>  - 


2u 


(37) 


The  result  for  the  collimated  beam  is  identical,  and  in  fact  could  be 
deduced  for  an  arbitrary  beam  focus  (Eq.  27)  through  conservation  of  energy: 


<I(£)>  = (^l)  J dL<|u(r)|‘ 
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\ 


(38) 


Thus  the  mean  irradiance  at  the  receiver  (illuminator)  plane  is  uniform  and 
independent  of  turbulence  level. 

B.  Covariance  and  Variance  of  Irradiance  at  Receiver 

The  correlation  function  of  the  irradiance  at  the  receiver  is 

Bi(£1,£2)  - <I(£1)I(£2)>  " <U(£1)U*(£1)U(£2)U*(£2)>  (39 


We  now  point  out  that  the  statistics  of  the  incoherent  reflection  from  a 
diffuse  object  are  gaussian,  from  the  Central  Limit  Theorem,  with  a phase 
which' is  uncorrelated  with  amplitude  and  uniformly  distributed  from  0 to  2tt. 
In  the  near  field  of  the  target  illumination,  we  may  expect  the  gaussian 
nature  of  this  field  to  apply  after  perturbation  by  turbulence.  A determina- 
tion of  the  exact  conditions  necessary  for  this  assumption  to  hold,  and  the 
generalization  of  the  solution  to  cases  where  it  does  not,  will  be  reported 
in  future  work.  The  gaussian  assumption  yields 


B = <U(£1)U*(£1)XU(£2)U*(£2)>  + <U(£1)U*(£2)XU*(£1)U(£2)> 


- <I(£1) ><i(£2)>  + 


l r (£x  .£2)  1 2 


(40) 
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where  T is  the  mutual  coherence  function  at  the  receiver.  It  follows 
that  the  covariance  of  intensity  given  by 


ci<£1»£2)  " bi<£1»£2^  ~ <I<£1) ><i(£2)>  " lr<£1*£2)|2 


(41) 


The  mutual  coherence  function  at  the  receiver  can  be  described  from 

13 

the  extended  Huygens-Fresnel  formula. 


2 

* (Sl)  Jj^  < U<£l>UV  > j »[*!<£!•%> 


- r2(£2»£2)]  { < exP  [’l'2(£1*£1)  + ’l'2*(£2»£2)]>  (42) 


where  R2(£2.£2)  are  the  distances  from  ^ to  ^ and  J>2  to  £, 

respectively. 

By  the  Fresnel  approximation 


Rl^-l’£l^  “ R2^-2,£2^  = 2L 


2 2.2  2 

* Pl_P2+Pl“P2  £1--1  £2*£2 


(43) 


Finally,  from  equations  (43)  and  (44): 
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Since  the  wave  is  incoherent  after  refl  ction  from  the  diffuse  target, 
the  coherence  function  at  that  plane  can  be  represented  by  the  Dirac  delta 
function: 


<U(£l)U*(p2)>  * <I(P1)>6(£1-£2) 


(45) 


Using  this  in  Eq.  (44),  r(£^,£  ) can  be  simplified: 


i’(£1»£2)  ' ( 2ttL  ) 


exp  [ 
5/3 


ik(pj-p2) 

2L 


]j"d£  <I(£)>  exp  y-  (£1“£2)*£| 


po' 


(46) 


In  the  absence  of  turbulence,  this  equation  is  entirely  identical  to  the 

14 

Van  Cittert-Zernike  theorem  of  coherence  theory,  which  is  identical  to  a 
result  obtained  by  Goodman  for  the  mutual  coherence  function  of  a pulsed 
optical  radar.  ^ 

To  complete  the  solution,  we  utilize  the  mean  intensity  at  the  target 
for- the  focussed  case,  Eqs.  (34)  and  (36): 


o 


We  thus  have: 


14.  Born,  M,  and  Wolf,  E.,  Principles  of  Optics.  Pergamon  Press,  N.Y.,  1975. 

15.  Goodman,  J.W. , Proc.  IEEE  5^,  1688-1700,  Nov.  , 1965. 
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Finally,  the  normalized  covariance  function  of  irradiance  for  the 
focussed  case  can  thus  be  written: 


(51) 


and  the  normalized  variance  is 


aT2  = 1.  (52) 

f 

For  the  collimated  case,  the  same  variance  is  obtained,  and  the  covariance 
is 


(53) 


In  general,  the  normalized  irradiance  variance  for  an  incoherent  wave  is 
known  to  be  unity.  The  probability  distribution  for  the  irradiance  is  the 
exponential  distribution,  and  for  the  amplitude  it  fits  the  Rayleigh  dis- 
tribution. 

The  covariance  scale  lengths,  for  the  focussed  case  (Eq.  51),  are 
and  pq,  although  the  present  analysis  is  clearly  applicable  only  when  the 
aQ  term  dominates  (po>ao>.  For  the  collimated  case,  the  dominating  scale  is 

^ (<ao).  These  terms  represent  the  "speckle"  sizes  for  focussed  and 

o 

collimated  beams  respectively. 

These  results  may  be  readily  extended  to  power  spectra  and  receiver 
aperture-smoothing  analyses,  which  will  be  given  in  a future  report. 
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